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Abstract: We show how the Dijkgraaf-Vafa matrix model proposal can be ex- 
tended to describe five-dimensional gauge theories compactified on a circle to four 
dimensions. This involves solving a certain quantum mechanical matrix model. We 
do this for the lift of the M = 1* theory to five dimensions. We show that the 
resulting expression for the superpotential in the confining vacuum is identical with 
the elliptic superpotential approach based on Nekrasov's five-dimensional general- 
ization of Seiberg-Witten theory involving the relativistic elliptic Calogero-Moser, or 
Ruijsenaars- Schneider, integrable system. 
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1. Introduction 



The Dijkgraaf-Vafa approach for calculating the superpotential of a four- dimensional 
super symmetric gauge theory via a matrix model [1-3] is by now well established. 
The matrix model seems capture all the holomorphic information of an Af = 1 theory 
in a lower dimensional — in this case a zero dimensional — system. Recently Dijkgraaf 
and Vafa [4] have proposed a much more general form of their approach which would 
apply, for instance, to higher dimensional theories compactified to four dimensions. 
This far-reaching proposal is very striking and without doubt deserves detailed study 
and independent testing. 

From the original proposal involving matrix integrals, it has turned out that 
relevant perturbations of the Af = 4 theory, the so-called Af = 1* theory, [3, 5, 
6], and its Leigh- Strassler deformations [7,8], can be solved since the associated 
matrix models were studied in other contexts [9, 10]. For the Af = 1* theory the 
matrix model results agree with the elliptic superpotential approach which is an 
alternative approach for calculating the holomorphic condensates [13,14]. In order to 
investigate the recent higher- dimensional proposal, it is natural to consider relevant 
deformations of the M = 4* theory lifted to five dimensions where the additional 
dimension is compact. The recent work of Dijkgraaf and Vafa would have us consider 
a quantum mechanical matrix model rather than a matrix integral. We shall solve 
this model in Section 2 and hence compute the superpotential in the confining vacuum 
of the theory. The results in the other massive vacua follow by a straightforward 
generalization of [6]. 

Fortunately we have an independent and highly non-trivial check on the result 
based on the elliptic superpotential approach adopted in [13,14]. The generalization 
to the five dimensional theory has not been described previously but it not difficult 
to establish the form of the exact elliptic superpotential in this case, generalizing the 
one for the Af = 1* theory in [13]. The point is that the Seiberg-Witten theory of the 
compactified five-dimensional theory was considered by Nekrasov [15]. What happens 
is that the Calogero-Moser integrable system which underlies the four dimensional 
theory is replaced by its relativistic version: the so-called Ruijsenaars-Schneider 
system [16]. The exact elliptic superpotential for the simplest mass deformation of 
the Af = 2 theory is identified with the basic Hamiltonian of this integrable system. 
Hence we are able to independently calculate, for instance, the gluino condensate in 
the confining vacuum of the compactified five dimensional theory and the result is 
entirely consistent with the expression calculated via matrix quantum mechanics. 

As a side remark, we shall find that the quantum-mechanical matrix model can 
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be reduced to a matrix integral and this integral is very similar to the one that 
described the relevant deformations of the Leigh- Strassler theory in [7]. It is rather 
natural, therefore, that the elliptic superpotential side of the story should involve 
the same integrable system [8]; namely, that of Ruijsenaars and Schneider. Surely 
this is a strong hint of some deeper connection between the integrable system and 
the matrix model. 



2. The quantum-mechanical matrix model 

We begin by formulating the superpotential of the compactified five- dimensional 
theory as a quantum mechanical matrix model following Dijkgraaf and Vafa [4]. From 
the point-of-view of the four-dimensional theory there are 3 adjoint chiral fields $j, 
i = 1, 2, 3. In the five- dimensional theory, the imaginary component of, say, $3 is re- 
interpreted as the component of the gauge field along the extra compact dimension, 
while the real part of $3 is the real scalar of the five- dimensional theory. We will 
take the compact dimension to have length f3. The other fields $1 and $2 form an 
adjoint- valued hypermultiplet. The superpotential of the effective four-dimensional 
theory is determined by the quantum mechanical system involving the fields 
The partition function of the quantum mechanical system involves the functional 
integral: 

3 

Z= I f[m] exp ( - {fig,)- 1 I dtW[*i]) . (2.1) 
J i=i J 

The action of the matrix model is a generalization of the one that describes the 
M = 1* deformation of the four dimensional theory: 

W[$i] = Tr(i<$> 1 D<$> 2 + m$ 1 <S> 2 + {J l cosh(P<$> 3 )) (2.2) 

where the covariant derivative is D§ 2 — d t $ 2 + [^3, ^2]- I n fact we have been a bit 
implicit in writing down ( |2.1| ) because we have not specified the measure. Part of the 
Dijkgraaf- Vafa prescription involves interpreting the matrix integral in a holomor- 
phic way. So the complex fields $i(t) and $> 2 (t) are subject to a particular reality 
condition, namely <&{ = $2, or equivalently $1 + $ 2 and — $2) are Hermitian, 
from the point-of-view of the functional integral. In particular the measure for the 
latter combination of fields is the appropriate measure for Hermitian fields. The 
field $3(t) is treated in a somewhat different manner as described in [4]. First of all, 
local gauge transformations can be used to gauge away the non-constant part of the 
component of the gauge field along the compact direction. This leaves large gauge 
transformations which shift the eigenvalues of the gauge field by integer multiples 



-2- 



of 2ti/ (3. This means that the natural variable is not the gauge field, but rather its 
holonomy around the circle expi(3A t . In the holomorphic point-of-view the quantity 
iA t is then naturally complexified to $3 by including the real scalar and then the 
prescription of Dijkgraaf and Vafa [4] is to interpret the integral as being over the 
t-independent quantity 

U = exp /3$ 3 (2.3) 

thought of as a unitary matrix. Notice that the final term in ( |2.2|) is the natural 
generalization of the Tr $3 in the four-dimensional theory and incorporates the nec- 
essary periodicity $3 — > $3 + 2m/ /3 and will allow us to make contact with the 
generalization of the condensate u<i described in [15]. 

Now following the logic of [3,5], we integrate out the fields $i(t) and $2(2) since 
they appear Gaussian in (|2.1|). The complication is that we must now integrate out all 
the fourier modes of these fields. We end up with a pure unitary (zero-dimensional) 
matrix integral 

dJJ exp(-ff7VTr cosh(/3$ 3 )) 

det sinh \ (im + $3 ® 1 - 1 ® $3) ' 

where we used the identity 

det (ft + f) = det sinh Up<p)/(\p) ■ (2.5) 



Now we are ready to perform the large- N limit saddle-point evaluation of the 
remaining unitary integral around the critical point appropriate to the confining 
vacuum. As usual in a unitarity matrix integral we can diagonalize U and work in 
terms of its eigenvalues 

U ~ (e*V..,e*") (2.6) 

at the expense of introducing the unitary matrix version of the Vandermonde deter- 
minant: 

* = / n «. n sm ,:g- gfa ^ < - e «* *>) ■ (") 

Apart from the potential, this is precisely the matrix integral that appears in the 
solution of the six vertex model on a random lattice [10] and has already been 
employed to describe relevant perturbations of Leigh-Strassler deformations of M = 4 
theories. However, it is almost identical to the matrix model that was considered 
in [11] which arose from taking a dimensional reduction of M = 1 supersymmetric 
Yang-Mills in four dimensions to one dimension. The only difference is that the 
periodicity in the eigenvalues in that reference is along the real axis, however, this 
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is only a superfical difference because our integral is interpreted in a holomorphic 
way. In addition, our model is also identical to the different relevant deformation of 
the Leigh- Strassler theory considered in [12]. In order to solve the model, we will 
largely follow the approach of [7] (or [12]) which is tailored towards the Dijkgraaf-Vafa 
application. 

In the large- N limit, the eigenvalues (f>i form a continuum and condense onto cuts 
in the complex z-plane (actually the cylinder due to the identification z ~ z + 2m). 
One can think of these cuts as arising from a quantum smearing-out of the classical 
eigenvalues. For the confining vacuum all the classical eigenvalues are degenerate 
4>i = and so we expect a solution in the matrix model involving a single cut which 
we take to extend from —a to a. Notice that at this point we diverge from the Leigh- 
Strassler case where the cut is not symmetrical about z = [7]. The extent of the cut 
and the matrix model density of eigenvalues p(z) will be determined self-consistently 
from the saddle-point equation in terms of the 't Hooft coupling of the matrix model 
S = g s N. The saddle-point equation is most conveniently formulated after defining 
the resolvent function 

u(z) = \ , j\{<P)d<P = l. (2.8) 

This function is analytic in z and its only singularity is along a branch cut extending 
between [—a, a}. The matrix model spectral density p(<p) is equal to the discontinuity 
across the cut 

u{4> + * e ) — w(4> — ie) = — 27rip(0) , e [—a, a] . (2.9) 

In this, and following equations, e is an infinitesimal. The saddle-point equation 
expresses the condition of zero force on a test eigenvalue in the presence of the 
large- N distribution of eigenvalues along the cut: 

p sinh (f) 

— = uj((f>+ie) — ie) — u>(4> + i(3m) — — i/3m) , G [— a, a] . (2.10) 

This equation can be re-written in terms of the useful function 

G(z) = — * coshz + iS(u(z + £2) - U (z - if*)) . (2.11) 

From this definition, one may easily deduce that G(z) is an analytic function on the 
cylinder — it < z < it, where Imz = n and Imz = — 7r are identified, with two cuts 
[-a + a + ^] and [-a - a- ^]. This is illustrated in Figure 1. Note that 
everything is periodic in (3m — > f3m + 2tt so we can always choose — ir < f3m < n. 



In terms of G(z), the matrix model saddle-point equation Q2.10D is 



G(0+^P±ze) = G{<P-i^Tie) , <f>e[-a,a}. (2.12) 
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Figure 1: The region over which the function G(z) is defined with its two cuts. The lines 
Im z = ±7r are identified. 

This equation can be viewed as a condition which glues the top (bottom) of the 
upper cut to the bottom (top) of the lower cut thereby defining a torus with two 
marked points corresponding to the infinities Re z — ► ±00. The function G(z) is then 
uniquely specified by gluing condition Eq. (|2.12| ) and asymptotic behaviour at large 
±Rez 

lim G(z) -> ^— e ±z + Oie TZ ) , (2.13) 

Re*-±oo V ' 4sin(^) V 7 ' ^ J 

which is a consequence of Eq. ( |2.11| ). 

The auxiliary torus E f is specified by a complex structure f which can be thought 
of as a function of the parameter a specifying the length of the cut. We can uniformize 
the torus by establishing a map to the complex w-plane (u being an auxiliary variable) 
quotiented by a lattice with periods 2iri and 2mr. The torus is shown in Figure 2, 
where for illustrative purposes f has been taken to be purely imaginary. As shown 
in Figures 1 and 2, the contour A enclosing the cut [— a + ^p, a + ^4p] anti-clockwise 
maps to one of the cycles of the torus while the contour B joining the two cuts maps 
to the conjugate cycle. 

The map z(u) from the w-plane to the z-plane is specified by the requirements 
that going around the contour A returns z to its original value, while traversing the 
contour B causes z to jump by an amount i/3m, which is the distance between the 
two cuts in the z-plane. Both these operations leave G unchanged implying that it 
is an elliptic function on the u-plane. Thus 

A : z(u + 2ni) = z(u) ; G{z(u + 2m)) = G{z(u)) , (2.14a) 

B : z(u + 2nif ) = z(u) + i(3m ; G(z(u + 2m?)) = G(z{u)) . (2.14b) 
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Figure 2: Uniformizing map from the z-plane to the torus Ef. 
This determines the following unique map z(u) from the w-plane to the z-plane: 

expz(u) = -f-^ . (2.15) 

0i(f + ^|f) 

The only singularities of G(z(u)) are simple poles at u — ±^rp corresponding to 
large ±Rez in order to incorporate the asymptotic behaviour ( [2.13| ). The fact that 
G(z(u)) is elliptic, along with the singularity structure, specifies it uniquely: 



G(z(u)) 



2sin^ 



e'Mf) 



(C(u - f ) - C(« + f ) + 2C(z/3) - f C(7ri)) . (2.16) 



In ( p . 1 5|) and (|2.16|) , the quasi-elliptic function ((u) is defined on the torus E T and 
8i(x\t) is a Jacobi theta function (see [17] for definitions and Appendix A for some 
useful identities). As expected from ( |2.14b|) , G(z{u)) is an elliptic function of u 
with two simple poles in the w-plane; z{u) on the other hand is only quasi-elliptic. 
Having determined G(z) in the elliptic parameterization we can now implement the 
Dijkgraaf-Vafa proposal in order to compute the superpotential in the confining vac- 
uum. 



According the Dijkgraaf-Vafa proposal, the gluino condensate of the gauge theory 
gets identified with the 't Hooft coupling S of the matrix model. From Eqs. (|2.S| ) 
and ( |2.11|) , the integral of G(z) around the contour A is equal to —2irg s N = —2ttS. 
Under the map to the torus this becomes an integral around the A-cycle: 



2mS 



/ / G(z(u))^P-du . 
du 



(2.17) 



The second ingredient required to determine the QFT superpotential is the variation 
of the genus zero free energy T§ of the matrix model in transporting a test eigenvalue 
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from infinity to the endpoint of the cut. This is obtained by integrating the force on 
a test eigenvalue, which can be expressed in terms of the function G(z) as 



-i(G{z + G(z , (2.18) 



from infinity to the original cut [—a, a}. This can be written as an integral over G(z) 
alone along a contour starting at the lower cut, going off to infinity and then back 
to the upper cut. This can be deformed into the contour B running from the lower 
cut to the upper cut as in Figure 1. Under the map to the torus this becomes an 
integral over the .B-cycle: 

m -if G{z{u)) d 4^du . (2.19) 



OS J B du 

The effective superpotential in the confining vacuum is obtained by extremizing the 
following expression with respect to S: 

dJ- 

W eS = N-^- - 2mrS , (2.20) 

where r is the bare coupling of the theory not to be confused with the complex 
structure of E?\ we shall shortly relate the two. 



Both integrals (|2.17|) and ( p.!9| ) are evaluated using standard elliptic function 
identities: 

^f^-Hf), (2,1) 



where 



It follows that 



so that 



siri /|i 0J( O |f) 



dW cS n . r 

T 



dS N 



(2.23) 



As a test of our expression for the superpotential we can consider the limit (3 — > 
and relate it to the known result [3,5,13]. In this limit, from ( 2.24J) we find 



WeS - -Np - ^fE 2 {^)(3 2 + ■ • • , (2.25) 



where E^ir) is the second Eisenstein series. This result agrees exactly with the 
expected result [3,5,13] since in this limit the potential behaves as V r ($) — > fi(l + 

k? 2 $ 2 + ■■•)■ 
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It is possible to use the quantum mechanical matrix model to calculate the 
superpotential in all the massive vacua of the theory. The idea is to consider special 
multi-cut saddle-point solutions as described in [6]. We do not include the details 
here, but they can straightforwardly be deduced from the aforementioned reference. 
The final result for the p th massive vacuum, where the confining corresponds to p — 1 
and the Higgs to p — N, is 



Recall that the result for the superpotential of the matrix model calculation in the 
four-dimensional M = 1* case does not have exact modular symmetry. For instance, 
the S transformation r — > — 1/r only relates the superpotential in the Higgs and 
confining vacua up to an additive anomaly. In the five-dimensional case described 
above we see that there is also a modular anomaly in the matrix model result but 
now of a multiplicative form. Interestingly under a general modular transformation 

ar + b 



ct + d 

if3m transformations like a point on the torus with complex structure r: 



(2.27) 



ifim . (2.28) 

cr + d v ' 



3. The Elliptic Superpotential 

Another method that has been used to calculate the exact values of the superpotential 
in all the massive vacua of the four- dimensional M = 1* theory is to compactify it 
on a circle of finite radius [13,14]. The effective superpotential is then a function 
of the dual photons and Wilson lines of the abelian subgroup U^l)^ -1 C SU(iV). 
These comprise N — 1 complex scalar fields X a , a = 1, . . . , N (with X a = 0) 

which naturally live on a torus of complex structure r because of the periodicity 
of each dual photon and Wilson line. The superpotential describing the M = 1* 
deformation is therefore constrained to be an elliptic function of the complex scalars 
X a . It turns out that this superpotential is identified as the basic Hamiltonian of 
the elliptic Calogero-Moser integrable system where the X a are position coordinates. 
The resulting superpotential, as originally found in [13] can be 

Elliptic = /i( \ P l - ^ E ~ Xfe )) > t 3 ' 1 ) 

where p{z) is the Weierstrass function defined on a torus with periods 2m and 2nir. 
Notice that the momenta can trivially be integrated out. What is particularly useful 
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about this superpotential is that it is independent of the compactification radius, 
and, therefore, yields results that are valid in the four-dimensional limit. In addition 
unlike the matrix model approach it encodes all the vacua of the M = 1* theory, 
both massive and massless, in a single superpotential. 

The question is how this elliptic superpotential is generalized in the five- dimensional 
theory compactified on a circle. In order to motivate the answer we need to recall in 
more detail why the elliptic Calogero-Moser integrable system appears in the basic 
M = 1* case. The reason is that, if we think of the A/" = 1* theory in terms of a de- 
formation by the mass term fi of the M = 2* theory, then the Coulomb branch of this 
latter theory is described by the moduli space of a Seiberg-Witten curve. This curve 
is precisely the spectral curve of the elliptic Calogero-Moser system [18]. Now for 
the five-dimensional theory, Nekrasov has made a conjecture for the Seiberg-Witten 
curve: it is precisely the spectral curve of a one parameter deformation of the elliptic 
Calogero-Moser system known as the elliptic Ruijsenaars-Schneider system. This is 
sometimes known as the relativistic elliptic Calogero-Moser system. 

Following the approach of [13, 14] it is natural to conjecture that the elliptic 
superpotential of the five-dimensional theory, for the simplest deformation fi cosh 
is the basic Hamiltonian of the Ruijsenaars-Schneider system. Based on this, and 
with reference to [15], we expect 



for some constant C. It is easy to see that the critical points which describe the 
confining vacuum of the M = 1* case 



massive vacua described in [13]. 

We can now compare the result of the quantum mechanical matrix model fl2.24j ) 
with the result from the elliptic superpotential: 

N-l 



Amazingly the results agree by virtue of the elliptic function identity (|A.8|) proved 
in the Appendix 




(3.2) 





(3.4) 



a=l 





N 



N-l 



(»)"->( 



) 



J] yj p (i/3m\T) - p(^\r) , (3.5) 



o=l 
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where we have emphasized the complex structure associated to the elliptic functions 
on the right-hand side. 



Appendix A: Some Properties of Elliptic Functions 



We use the (quasi-) elliptic functions p(u), ({u), Oi(%-\t) associated to a torus of 
periods 2a>i = 2ni and 2uj 2 = 2i\ir as defined in [17]. An important equation relating 
them is 

We also use the heat equation 

d 2 9 i (x\r) 4 dOAxW) /A s 

'V - + - V =0 A.2 

ox 1 1x1 or 

and the relations 



7T 



2 



vr 2 0i"(O|r) 



CMu! = — ■ e 2 (t) = n \ ; / . (A.3) 

sv u 1 12 y ! 12 01(O|r) V ; 

In the remainder of the Appendix, we establish a crucial identity. We start 
with the following identity established in [7] (using ( |A.1| ), ( |A.3| ) and other standard 
properties of elliptic functions) 

a r u m ~ a r u i^ + r 2 ^ / l^v tv t~ C( jv 'J • ( A - 4 ) 



We now integrating this expression, using the relation (|A.1Q on the right-hand side, 
to arrive at 



9 ' " ' T 1 v ~ : 



• 2* 1 JV> 



nA(f-iw(f ( a - 5 ) 



Further, on the right-hand side, we can employ the relation 







D 


r) 






T 


? 



p(v) - p{u) (A.6) 



to derive the identity 



N-1 



^i(ii^) = ^(akr n wwypM - pc 23 ^) • ( a - 7 ) 

0=1 

From the ratio of this, and its derivative evaluated at u = 0, we have our crucial 
identity 
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